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A Global Study of the Behaviour of Black Hole X-ray Binary Discs 
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ABSTRACT 

We investigate the behaviour of the accretion discs in the outbursts of the low-mass black- 
hole X-ray binaries (BHXRB), an overview of which we have presented previously. Almost 
all of the systems in which there are sufficient observations in the most disc dominated states 
show a variation of the disc luminosity with temperature close to L oc T 4 . This in turn implies 
that in these states, the disc radius, R[ n , and the colour correction factor, f co \, are almost 
constant. Deviations away from the T 4 law are observed at the beginning and end of the most 
disc dominated states, during the intermediate states. Although these could be explained by 
an inward motion of the accretion disc, they are more likely to be the result of an increase 
in the value of / co i as the disc fraction decreases. By comparing the expected and observed 
disc luminosities, we place approximate limits on the allowed distances and masses of the 
BHXRB system. In a number of cases, the measured distances and masses of the BHXRB 
system indicate that it is possible that the black hole may be spinning. 



Key words: 

binaries 



accretion, accretion discs - binaries: general - ISM: jets and outflows - X-rays: 



1 INTRODUCTION 

The outbursts of black hole X-ray binaries (BHXRBs) are dra- 
matic and intriguing events. They have the potential for allow- 
ing the study of the physical and emission processes close to the 
event horizon. The accretion process and associated intermittent 
jet-production results in emission across the electromagnetic spec- 
trum. In the study presented here, we focus on the X-rays, as these 
arise from the inner parts of the accretion disc and flow. At the other 
end of the spectrum, the radio emission is thought to arise from a 
synchrotron emitting jet. Therefore the radio emission is a good 
tracer of whether a jet is active or not, and the X-rays are good at 
determining the state of the accretion flow. 

In the now commonly-accepted picture of the changes that oc- 
cur in the BHXRB system, the BHXRB spends most of its time 
in a quiescent state. There the total luminosity of the BHXRB is 
very low, in all bands. As the outburst starts, the X-rays are char- 
acterised by a hard emission spectrum with a powerlaw slope of 
r ~ 1.5 - the "hard state". As the X-ray luminosity r ises, the ra- 
dio luminosity rises in step JCorbel et alj|200d, l2003h . The radio 
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spectrum also indicates the presence of a steady jet emitting syn- 
chrotron radiation. As the outburst progresses, the disc spectrum 
becomes increasingly dominant, eventually softening the entire X- 
ray spectrum as the BHXRB enters the "soft-state". This transi- 
tion is very fast compared to the speed at which the luminosity 
rose. Over the course o f weeks or even months, t he disc luminosity 
and temperature decay l Gierlinski & D one 2004), as the disc dom- 
inance decreases founn et alj|2010h . Eventually the source returns 
to the "hard-state" and the luminosity continues to fade. For further 
details on this pictu r e of the progression of BHXRB outbursts see 
iFender et alj d2004h:li5one & Gierlinskil d2003h;lHoman & Belloni 
d2005h: iRemillard & McClintockl feOOfJ) ; iDone et all d200' 
iBellonJ d2010h . 
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The accretion onto the compact object is what drives their lu- 
minosit y in X-rays, which a rise from the disc, the corona and even 
the jet ( iRussell et al.ll201oT) . Studying emission from the disc al- 
lows the accretion process, and also the behaviour of the mate- 
rial within the disc as it approaches close to the compact object 
to be investigated. The st andard theory of accretion discs from 
Shak ura & Sunvaevl Jl973h shows that the accreting material will 
form a geometrically thin, but optically thick disc, the inner extent 
of which depends on the spin of the black hole. 

Previous studies on the behaviours of the BHXRB accre- 
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tion discs have selected those observations where the disc was 
dominant, i n order to ensure that the disc parameters were well 
determined dGierlifiski & Donell2004 bone et al.ll2007h . This al- 
lowed a detailed study of the disc emission from well characterised 
BHXRBs to be carried out. Using the full archival coverage of the 
RXT^\ satellite, which has been observing BHXRBs for 13 years, 
we present a study o f all observa t ions i n which a disc was detected 
using the analysis of iDunn et alj d2010l) . This allows us to investi- 
gate the properties of the disc in the transition periods, between the 
fully disc-dominated and powerlaw-dominated states, as well as in 
the disc dominated states. 

In Sections [2] and [3] we r ecap t he data reduction proceedure 
presented in lDunn et alN2008l . l2rJlCf) and the final BHXRBs which 
were selected for this study. The behaviour of the disc's temper- 
ature and luminosity are discussed in Section [5] The deviations 
from the expected behaviour of the disc luminosity and tempera- 
ture are presented in Section [6] where the inner radius of the disc 
is investigated, and Section [7] where larger departures are linked 
to the colour temperature correction. In Section [8] we investigate 
the limits which can be placed on the distances, masses and spins 
of the BHXRBs from the observations. The degeneracy of the bro- 
ken powerlaw model with the disc model, and the evolution of the 
powerlaw in the hard state are presented in Sections [9l andl 1 01 



2 DATA REDUCTION 

In the analysis of the disc properties of the sample of BHXRBs 
presente d in this work we use the analysis of RXTE data de- 
ta iled in IDunn et al. (2010). We recap the main points, but refer 
to lDunn et al.1 J2008L [20100 for more details. 

We use all the available data publicly available in the RXTE 
archivfl This gave a baseline of around 13 years to study the evolu- 
tion of the disc properties during the numerous outbursts observed 
within that time. All data were subjected to the same data reduction 
procedure, in order to minimise differences arising from different 
data reduction routines. 

Both the Proportional Counter Array (PCA) and High Energy 
X-ray Timing Experiment (HEXTE) data were required when fit- 
ting the spectra, as the HEXTE data allows the powerlaw to be con- 
strained at high energy when the PCA data are dominated by the 
disc. We followed the procedure outlined in the RXTE Cookbook 
using the tools from HEASOFlQ version o.o.z. 

To reduce variations the between observations further, we only 
use the data from Proportional Counter Unit (PCU) 2 on the PCA 
as this has been on throughout the RXTE mission. Our analysis con- 
centrates on the bright periods when the BHXRBs are in outburst, 
and so we use the bright model background for all data. Lower 
count rates are more likely in the inter-outburst periods, and so this 
choice of a single background is unlikely to bias our results. 

In order to proceed with the spectral fitting, we require a PCA 
observation with at least 1000 background subtracted counts, and 
a HEXTE observation with either Cluster A or B (or both) with at 
least 2000 background subtracted counts. The other HEXTE cluster 
has to have at least a positive number of count jf| This count restric- 
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tion is in place to try to ensure that the spectra which are fitted are of 
good quality and fit within a reasonable time with well-constrained 
parameters. 

The spectra were fitted in XSPEC (vl2.5.0an). In order to study 
the disc parameters in detail, we needed to analyse the spectra to the 
lowest energies possible. The relation between the channel num- 
bers of the PCA instrument and the energies they correspond to 
has drifted over the 13 years of the mission. However, all channels 
below number 7 are not well calibrated for spectral analysis. We 
therefore choose to ignore PCA channels ^ 6, which corresponds 
to around 3 keV, but the exact energy has drifted over time (see the 
RXTE documentation). We also ignore PCA data > 25 keV, and 
HEXTE data < 25 keV and > 250 keV. 

To be able to characterise the state of the BHXRBs as they 
go through an outburst we fit three types of base model - unbroken 
powerlaw (POWER, PL), broken powerlaw (BKNPOWER, BPL) and 
powerlaw + disc (POWER + DISKBB, DPL). These allow the study 
of the non-thermal component using the POWER/BKNPOWER pa- 
rameters, and the disc using the DISKBB parameters. To study the 
presence and change in the iron line we add an optional 6.4 keV 
gaussian feature to all these spectra, giving in total six models 
which were fitted. The low energy sensitivity of RXTE is insuffi- 
cient to allow the Nh to be determined from the spectra, and so we 
fix this value to the accepted value for each BHXRB (see TableQ}. 

From the six fitted models, we select the best fitting one on 
X 2 terms. However if this is not the simplest model, we then deter- 
mine whether the increase in complexity of the model is significant 
using an F-test with V < 0.001 as th e significance level. For the 
complete routine see lDunn et alj ( l2010h . but a quick outline is de- 
scribed below. When the best fitting model is complex but contains 
no gaussian component, we test this best fitting model against the 
simple powerlaw result. If the best fitting model is complex and 
contains a gaussian component we first test whether the underlying 
complex continuum model is an improvement over the simple pow- 
erlaw, and if it is we test whether a line is required in this complex 
model. When the complex continuum is not an improvement over 
the simple powerlaw a num ber of further steps are performed, as 
detailed in lDunn et al] d2010h . 

Once the best fitting model has been selected, we further cut 
the observation number by removing any observation whose 3 — 
10 keV flux is less than 1 x 10 -11 erg s _1 , where the flux was not 
well determined or where the powerlaw was not well constrained 
(even if the disc was). The flux cut was performed to focus on the 
periods in which the BHXRBs are in outburst, and so streamline the 
data reduction process. We also removed those fits whose \ 2 > 5.0 
as these are spectra which are not well fit by any of the models 
available within our automated procedure. The distribution of the 
X 2 of the best fitting models is shown in lDunn et alj ( 12010 ) Fig. 
2. The majority of fits are clustered around x 2 ~ 1, but there is a 
large tail to higher values. As the spectral fitting in this work has 
been automated, such large tail is expected. 

2.1 Model fitting issues 

The relatively high lower energy bound for the RXTE response lim- 
its our ability to detect discs when they are not dominant. The maxi- 
mum power emitted by the DISKBB model occurs around 2.4fcTbi S c 
which is usually around the lower limit of the PCA bandpass (for 
d^scj; at ~ 1 keV). Therefore we rarely detect the peak of the disc 
emission, and more usually observe the Wien tail. Using a sim- 
ple powerlaw to model the non-thermal continuum, even when in- 
cluding HEXTE data, does not allow for small breaks or curvature 
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within this component. If a disc component was included in these 
observations, it was found to try and fit these small curvatures in the 
powerlaw rather than any true underlying disc component, result- 
ing in unphysical disc parameters. We therefore limited the mini- 
mum temperature for the disc during the fitting to feT = 0.1 keV. 
Furthermore we then penalise the \ 2 of any model which has a 
fcaT < 0.4 keV when selecting the best fitting model. We note 
that in doing this we are limiting our sensitivity to low tempera- 
ture discs, in the intermediate and hard states for example, and are 
probably excluding a few accurate disc fits. We investigate further 
degeneracies between the disc and broken powerlaw models in Sec- 
tion|9] 

More complex models, for example Comptonization, would in 
principle give more information on the state of the system in these 
disc dominated states, as it links the non-thermal emission to the 
disc temperature. However, in order to freely fit all the parameters 
of the Comptonization models a high signal-to-noise observation 
is required. Not all of our observations have sufficient counts to be 
able to do this; in fact very few would allow all parameters to be 
determined from the observations. Although fixing some parame- 
ters would allow these models to fit successfully, this goes against 
the methodology of this work, by a priori constraining parameters 
differently for different states. 



3 SELECTED OBJECTS 

After all the data reduction, dead-time and selection the 15 Ms of 
raw RXTE data was trimmed to ~ 10 Ms in 3919 observations, with 
well fitted spectra and high enough fluxes and counts. The sample 
of objects was not designed to be complete in any way. We selected 
objects which were well known BHXRBs in the literature as well 
as those which were known to have outbursts which had been well 
monitored by RXTE. The set of BHXRBs analysed in this sample, 
along with their physical parameters (where known) and the final 
number of observations used in this study are shown in Tables Q] 
and|2] 

There are two notable BHXRBs which were purposely not in- 
cluded in this study (e.g.Cyg XI and GRS 1915-105). These two 
sources were not included for a number of reasons. One was a 
purely practical one resulting from the shear amount of data avail- 
able for these sources. The reduction of all the observations in the 
scheme outlined above would h ave dominated a ny of the global 
studies presented both here and in lDunn et al.l d2010l) . and selecting 
certain parts would have gone against the philosophy of the study, 
by not including all of the available data. Secondly, the behaviour 
of these sources is not easily explained by the outburst model pre- 
sented in iFenderetal J ( 120041) . In the following sections, we use this 
outburst scheme and the states it describes to explain the behaviour 
of the disc and powerlaw components. As the behaviour of these 
two well studied BHXRBs do not easily fall fit into this scheme, 
we actively decided to not include them in the study at this time. 

Many of the masses and distances are unknown or not very 
well constrained. Where they are unknown we have assumed values 
of IQMq and 5 kpc respectively. These uncertainties effect the cal- 
culation of the Eddington Luminosities (Z/Edd) for these BHXRBs 
which are used extensively throughout this analysis to scale the 
BHXRBs to one another. Until the distances and masses are well 
determined, there will always be some uncertainty when compar- 
ing between sources. 

The investigation presented in this work concentrates on the 
variation of the disc characteristics during the outburst as the 



Table 2. Observation Numbers, Times and Disc Detections 



Object Selected Obs Exposure Disc Detections 



Ms 


41 J 1 S43-47 


61 


0.147 


36 


4U 1630-47 


704 


1.371 


491 


4U 1957+115 


59 


0.260 


25 


GRO J1655-40 


484 


1.829 


368 


GRS 1737-31 


5 


0.045 


2 


GRS 1739-278 


g 


0.017 


6 


GRS 1758-258 


g 


0.007 


9 


GS 1354-644 


g 


0.049 


1 


GS 2023+338 





0.000 


o 


GX 339-4 


709 


1.682 


284 


H 1743-322 


346 


0.998 


224 


XTE Jl 1 18+480 


81 


0.170 


1 


XTEJ1550-564 


365 


0.833 


168 


XTE J 1650-500 


108 


0.191 


37 


XTEJ1720-318 


63 


0.125 


33 


XTEJ1748-288 


21 


0.074 


12 


XTEJ1755-324 


2 


0.006 


1 


XTEJ1817-330 


123 


0.329 


100 


XTE J 1859+226 


121 


0.292 


101 


XTE J2012+381 


15 


0.036 


15 


LMC X-l 


69 


0.349 


64 


LMC X-3 


471 


1.048 


173 


SAX 1711.6-3808 


13 


0.029 


5 


SAX 1819.3-2525 


48 


0.114 


2 


SLX 1746-331 


28 


0.091 


14 


Totals 


3919 


10.09 


2172 



changes in the disc parameters are the most prominent changes in 
the spectrum during a BHXRB outburst. Our results are therefore 
dominated by those objects which have had outbursts well moni- 
tored by RXTE. Roughly this "removes" all the BHXRBs from our 
study which have only had a few RXTE observations. Some objects 
which have had a comparatively large number of observations are 
not observed to undergo the canonical outburst structure outlined 
in SectionQ] These sources are less able to show what changes disc 
undergoes during a complete outburst, but are still useful for the 
hard/powerlaw dominated states. 



4 DISC PARAMETERS 

The discs around bl ack holes are thought to be optically thick and 
geometrically thin dShakura & Sunvaevlll973l) . The spectrum ex- 
pected from this kind of disc around a non-rotating black hole is 
easily calculated. It is the sum of a set of blackbody spectra, one 
for each radius, R, with a characteristic temperature T c s(R). The 
total spectrum resulting from this sum is then a multicolour disc 
blackbody, with a peak temperature T e g -, max coming from close to 
the innermost stable orbit. However, this spectrum is effected by 
the opacity of the disc, which results in a colour temperature cor- 
rection factor, / co i dShimura & Ta kahara 1995; Merloni et al. 2000; 
Davi s et al . 2006). This factor was shown bv lShimura & Ta kahara 
d 19951) tobe ~ 1.8 for almost all black hole masses and emis- 
sion luminosities and is discussed further in Section [7] Of course, 
the descrip tion of the disc may not be quite as simple as en- 
visaged by Sha kura & Sunvaevl d 1973b and radiat ively inefficient 
flows (e.g. Advection Dominated Accreti on Flows. | lchimaru! 197?] ; 
lNaravan & YiH l994) or slim discs (e.g. lAbramowicz et alj|l98a) 
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Table 1. X-ray Binary Parameters 



Object 




D 




Pcrb 


M* 


Inclination 




(Mq) 


(kpc) 


(xlO 22 cm- 2 ) 


(h) 


(Mq) 


(°) 



/ITT 1 C/11 /IT 

4U 1j43-4/ 


9.4 ± 2.0 


(1.2) 


7.5 i 0.5 


(3,4) 


0.43 


(2,4) 


26.8 


(4) 


2.45 


(1) 


21 


(2) 


/ITT 1 tin A 1 


M nl 
[10J 




1 n n _i_ k n 
1U.U ± O.U 


(5) 


> 6 


(o) 














/IT T 1 Q^"7_i_1 1 ^ 
4U iyj l+l 13 


[10J 








n 1 c 
u. lo 


( ') 


n q 


(oj 


i n 
l.U 








Ti tec /in 


7.0 i 0.2 


(10,11) 


3.2 i 0.2 


(4,12) 


0.8 


(13) 


62.9 


(4) 


2.35 


(10) 


70 


(48) 


/""D C 1 111 1 1 


M nl 
[10J 




Trr 1 
[ 5 ] 




p, n 
O.U 


(14) 














(jKIS 1 /3y-z/o 


ri nl 
[10J 




8.5 i 2.5 


(15) 


_' 


(15) 














(jKS 1 /jo-zjo 


ri nl 
[10J 




[5J 




1.50 


(16) 


18.5 


(17) 










c i i^/i t/i/i 
Ljj 13j4-044 


rr q inn I on 

> /.o — 1U.U ± Z.x) 


(1) 


> z ( — oo ± 


(18) 


O. ( Z 


(is, ty) 


Dl.l 


(IS) 


i no 
l.Uz 


( I ) 






/~* c inn i qio 
(jj ZUzi+iJo 


10 i 2 


(1) 


4.0 i 2.0 


fA\ 

(4) 


0.7 


(4) 


155.3 


(4) 


0.65 


(1) 






(ja i jy-4 


5.8 i 0.5 


(20) 


8.0 i 4.0 


(21) 


0.4 


(22) 


42. 1 


/ A \ 

(4) 


0.52 


(20) 


i n 
4U 




ri 1 /4J-JZZ 


M nl 
[10] 




Trr 1 
[ 5 ] 




O -1 

Z.4 


(23) 














Ytf; Ti 1 1Q j_ /ion 
A 1 E, J 1 1 lo+4oU 


a o 1 n A 


(1.24) 


1 *7 J_ n nc; 


(,ZJ,ZO) 


n m 
U.U1 


/o^~i 


a no 


(4) 


n oc 


1 I \ 


On 


(2b) 


TV IE J 1 J JU-JOt 


1U.U HI l.U 


\->) 


K Q _|_ O Q 

o.o m z.o 


(A\ 

W 


U.UO 


(Z/J 


O I . u 




1 ,; in 

l.OU 




79 
t z 


(•>) 


XTEJ1650-500 


< 7.3 = 6 ± 3 


(28) 


2.6 ±0.7 


(29) 


0.7 


(30) 


7.7 


(28) 






30 


(49) 


XTEJ1720-318 


[10] 


(31) 


> 8 = 8 ± 6 


(31) 


1.24 


(31) 














XTEJ1748-288 


[10] 




> 8 = 10 ±2 


(32) 


7.5 


(33) 














XTEJ1755-324 


[10] 




[5] 




0.37 


(34) 














XTEJ1817-330 


< 6 = 4± 2 


(35) 


> 1 = [10] 


(35) 


0.15 


(35) 














XTE J 1859+226 


10 ±5 


(36) 


6.3 ± 1.7 


(4) 


0.34 


(36) 


9.17 


(4) 


0.9 


(36) 






XTEJ2012+381 


[10] 




[5] 




1.3 


(37) 














LMC X-l 


10 ±5 


(38) 


52 ± 1.0 


(39) 


0.5 


(13) 


93.8 


(40) 






45 


(38) 


LMC X-3 


10 ±2 


(41) 


52 ±1.0 


(39) 


0.06 


(42) 


40.8 


(43) 


6 


(41) 


60 


(41) 


SAX 1711.6-3808 


[10] 




[5] 




2.8 


(44) 














SAX 1819.3-2525 


10 ±2 


(46) 


10 ±3 


(46) 


0.1 


(47) 


67.6 


(46) 






65 


(46) 


SLX 1746-331 


[10] 




[5] 




0.4 


(45) 















Many of the objects do not have well determined distances or masses. In this case we have taken the distances to be 5 kpc and the masses 10 Afm . A recent 

I ; — m i w 

critical look at the distance estimates for GRO J 1655-40 by Foellmi (2009) indicates a revised estimate of the distance of < 2.0 kpc. References: 

m lRitter & KoltU2003l) , (2) |Park et all J20041) , ( 3) IOrosz et aljj2002h (4)|jonker & NelemansU2004l) , (5) lAugusteiin et alj j200lh , (6) lTomsick et alj<2005l) 

. (7)lNowak et alj<2008l) ■ (8))Thorstensenl jl987h . (9)1 ShahbazeUujJ 19961) ■ (10)|Hvnes etai]h998h . (ll)IShahbaz et aljjl999l) (12)lHiellming &Ruperj 

il995l) . (13)lGierlihski et alJ j200lb ."(14) ICui et alj il997h . (15) iGreiner et alj ll996!) 7(16) IPottschmidt et alj <2006t) . (1 7) ISmithet alJ <2002l) (18 ) 

ICasares et aDi2004h (19)lKitamoto et alj jl990t) , (2 0) lHvnes et alj J2003I) , (2 1) lzdziarski et alj J2004l) . (22)lMiller et alj l2004h . (23)ICaDitanio et al J J2005h 

(24) Wagner et alj <200ll) (25)|Chatv et alj 120031) . (26)lGelino et alJ i200d) . (27)lGierlihski & DoneU2003l) ( 28)IOrosz et alj )2004h . (29)]Homan et alj 

J2006I) . (30)lMiniutti et aljj2004l) . (31)ICadolle Bel et alj l2004h (32)lHiellming et aljjl998l) (33)lKotani et alj j2000h . ( 34)lRevnivtsev et aljjl998l) ■ (3 5) 

ISala et alj J2007h (36)lHvnes et alj J2002I) (37)ICampana et 420020 , ( 38)lHutchings et alj f 19871) (39)ldi Benedettd Jl997l) (40)IOrosz et alj J2008l) 

(41)ICowlev et alj jl983l) , (42)lHaardt et alj j200ll) . (43)lHutchings et alj i2003h . (44)lin't Zand et alj J2002h . (45) IWilson et alj J2003h . (46) lQrosz et alj 

J200lh , (47) lin't Zand et al] j2000|) , (48) lvan der Hooft et alj ll998l) . (49) ISanchez-Fernandez et alj J2002l) 



m ay exist. However, we con centrate on the disc model proposed 

bv lShakura & Sunvaevl Jl973h in this st udy. 

We use the physic al description in iGierlinski & Pond J2004I) : 

iGierlinski et alj j 19991) to calculate the relation between the disc 
luminosity and the temperature as for a Schwarzschild black hole, 



Lois 



0.583 



1.8 

/col 



M 



kT, r 



IkeV 



(1) 



L Edd V /cot ) ViOMe 

which assumes a constant inner disc radius, R ln . We investigate the 
effects of the black hole spin in Section [8] We include the adjust- 
ments to the observed di sc temperature, Tobs, for relativistic effects 
close to the black hole l lGierlihski &~D one 2004). We add a 4 per 
cent temperature shift to acco unt for the str e ss-fre e boundary layer, 
and also the adjustment from IZhang et "all h997t) which accounts 
for the strong gravitational potential. 

Tmax = Tobs//GR(#, <!*)£., 

where £ = 1.04 is for the stress-free boundary layer, 8 is the incli- 
nation angle and a* the dimensionless spin parameter. Out of the 
~ 3900 observations ~ 2200 have disc detections (see Table O. 
Although the number of observations in which a disc is well de- 
termined depends on the state of the BHXRB at the time it was 



observed, we show the number of disc detections so that it is clear 
that a few BHXRBs have many more detections than most of the 
others. Therefore our conclusions are depend more on the results 
from these BHXRBs. 

In the data reduction routine and best fitting model selection 
proceedure we have been conservative in determining which ob- 
servations have discs (see Section |2~Tl >. Initially we investigate the 
degree to which the BHXRB discs follow the expected L — T rela- 
tion (EquationQJ, using an / co i = 1.8 and assuming that the inner 
radius of the disc is constant. In later Sections we relax these as- 
sumptions. 

In the following, we define the Powerlaw Fraction (PLF) and 
the Disc Fraction (DF) as 



PLF : 
DF = 



il-100kcV, PL 



Lo. 001-100 kcV, Disc + il-100kcV, PL 
-^0. 001-100 keV, Disc 



L 



001-100 kcV, Disc 



+ L 



1-100 kcV. PL 



following IPunnetalJ d20ld) as well as iDunn et alj d2008h : 
Kordi ng et aTl l2006). These two quantities, used when creating 
Disc Fraction Luminosity Diagrams, allow the natural separation of 
the outburst into two states - powerlaw and disc dominated. These 
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correspond roughly to the hard and soft states more commonly used 
in BHXRB studies. For a n in depth study of the relation between 
these state conventions see Dunn et 



4 



5 DISC TEMPERATURE AND LUMINOSITY 

We show the variation of the disc temperature with unabsorbed disc 
luminosity for each BHXRB individually in Appendix Fig. I A. 11 
The errorbars are only from the uncertainties arising in the spectral 
fitting. We do not include the uncertainties in physical parameters 
of the BHXRB system (e.g. mass and distance), as in many cases 
the physical parameters are unknown, and would further compli- 
cate the diagram. We also show the theoretically expected L — T 
relation for / co i = 1.8 for each BHXRB on each diagram in Ap- 
pendix Fig. IA. li as the dashed black line, using the masses as shown 
in Table [TJ (again without including the uncertainties). Also shown 
in Appendix Fig. I A. 1 1 is a schematic showing the motion of the 
BHXRB through the L — T plane as an outburst progresses. For 
clarity, for the remainder of this section the theoretical L — T re- 
lation is that from Eq. Q] under the assumption of a constant inner 
disc radius and colour correction factor. 

In Appendix Fig. I A. 1 1 it is clear that most of the BHXRB 's 
discs do closely follow the theoretically expected L — T relation. 
We fit the most disc dominated points (DF > 0.8) of each BHXRB 
with a powerlaw in the log 10 T — log 10 Z/Disc plane, where for a 
constant size black body a slope of four is expected. We show the 
resulting slopes in Table [3] and a histogram of their distribution in 
Fig- El The best-fit Gaussian distribution to the histogram peaks at 
a slope of 4.48. We select the most disc dominated observations 
in order to focus on those where the disc parameters (tempera- 
ture and disc normalisation in XSPEC) were very well determined 
and also to exclude points close to the intermediate state, where 
the relation may not a pply (similar to the selection performed in 
lGierliriskT& Done 2004). The behaviour of the disc temperature at 
smaller disc fractions, corresponding to states closer to the inter- 
mediate states, and are discussed in Section|7] 

Of the ten BHXRB s studied by IGierlinski & Pond ( 12004b . 
nine are included in our study. For four of the BHXRBs, their dis- 
tribution of observations match between their study and those in 
Appendix Fig [AT] (GRS 1739-278, XTE J2012+381, LMC X-l 
and LMC X-3). There are many more observations of GX 339- 
4 presented in our study, and so it is difficult to determine 
any differences between the two studies. Of the remaining four 
BHXRBs (GRO J1655-40, XTE J1550 -564, XTE J1650-500 and 
XTE J1859+226), the trends observed in lGierlinski & Pond ( 120041) 
show very clear and linear ~ T 4 relations. However, in Appendix 
Fig I A. 11 we find that, although the most disc dominated observa- 
tions do on the whole follow the expected L — T relation, there 
are a large number of points at low disc fractions which fall "be- 
low" the expected L — T relation (see also Section[7). In this study, 
we include and show all observations i n which a disc+powerlaw 
model was the best fit, whereas those in IGierlinski & Pond ( 12004) 
select "disc dominated spectra" where up to 15 per cent of the total 
bolometric emission can be present in a Comptonized tail. In Ap- 
pendix Fig lA.ll we only fit those observations for which PF> 0.8 
and therefore the apparent observed differences are large because 
of the plotted low disc fraction points. 

In a number of the BHXRBs, the statistically best fitting line 
is not similar to the expected L — T relation. This mismatch in 
between the slopes of the expected L — T relation and that of the 
best fit to the most disc dominated states could arise from the lim- 



2 3 4 5 6 7 

Exponent 

Figure 1. The distribution of the best fitting slopes in the most disc domi- 
nated states. The indicated Gaussian is a least-squares fit to the histogram, 
and peaks at 4.48 with a width of 0.33. GRS 1739-278, XTE J1650-500 
and LMC X-l are beyond the edges of the plot (see Table|5J- 



itations on the range of disc temperatures probed and the spectral 
response of RXTE. If only a few observations have a detected disc, 
then the variation in the disc temperature and luminosity may be 
small, which could mask a T 4 trend if the scatter is naturally high, 
for example GRS 1739-278 and LMC X-l. 

In a large number of cases, although the shape of the relation is 
close to that of the L — T relation, the normalisations are not always 
correct (see Equation[TJl. Given the uncertainties in the masses and 
distances an offset between the observed and expected behaviours 
is not unexpected. We also note that the expected L — T relation 
assumes that the black hole is not spinning. We discuss the effects 
of these unknown system parameters further in Section[8] 

Therefore it is clear that under the assumption of a constant 
/coi and a constant inner disc radius, R[ n , the majority of the 
BHXRBs show a variation of the disc luminosity, Lnisc 6c T 4 in 
the most disc dominated states over an order of magnitude in lu- 
minosity and a factor of tw o in disc temperature. This has been 
found by earlier studies (e.g. Gierlihski & Pone 2004; |Pavis et al.l 
2006; Pun n et alj|2008l) . though in some cases the analysis as re- 
stricted to the most disc dominated (soft) states. Specifically, the 
ratio of Ri n : f co \ is constant in these states. However, it would 
be very strange if the dynamics of the accretion disc and the radi- 
ation transfer processes conspired to keep this ratio constant, and 
so it is likely that each quantity is itself, constant, in the disc dom- 
inated state. However, what of the behaviour in the less disc domi- 
nated states? As noted above, many of the observations at low disc 
fractions fall "below" the expected "L-T" relation. We investigated 
whether this could result from slightly different L — T relations 
from different outbursts of the same BHXRB, however no clear 
or variation was found. A number of the low disc fraction observa- 
tions are seen as almost perpendicular deviations from the expected 
L — T relation. The trend for these "spurs" is for the disc luminosity 
to decrease as the disc temperature increases and are observed both 
at the beginning and the end of the outbursts (see e.g. 4U 1543-47 
and GRO J 1655-40). We discuss these deviations further in Section 

13 

There are a few cases where the best fit line to the most disc 
dominated observations is radically different to the expected L — T 
relation even though there are a large number of observations, e.g. 
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XTE J 1 650-500, LMC X- 1 . It is almost as if whatever is causing the 
"spurs" dominates the variation of the disc temperature and lumi- 
nosity in these objects. In these two BHXRBs there are two effects 
which conspire to give best-fit relations different to the expected 
L oc T 4 . Firstly, in both of these sources the majority (if not all) 
of observations have the same luminosity. In XTE J1650-500 those 
observations occur just after the disc fraction reaches 0.8, rather 
than in the most disc dominated state of the outburst (DF> 0.9). If 
we select those observations with DF> 0.9 as opposed to DF> 0.8 
then the slope of the best fitting line is 9.52. Whereas in LMC X-l 
the only a few observations have a sufficiently high disc fraction to 
be fitted, with almost no variability in the total luminosit)[fl Sec- 
ondly, as seen in a number of the other BHXRBs, observations at 
intermediate states, or those which do not have overly strong disc 
fractions show "spurs" running approximately perpendicular to the 
L — T relation. Both of these effects are combined with the scatter 
observed in the the relations exhibited by other BHXRBs. 

Therefore, if there are few points in the most disc dominated 
states (XTE J 1650-500) or there is a small variation in the lumi- 
nosity (LMC X-l), it is possible that the best fit relation will not be 
close to the theoretically expected one. However, in XTE 11650- 
500, the best fit relation to the observations with DF> 0.9 is much 
closer to the expected L — T relat ion. In a study of the discs in 
BHXRBs iGierlifiski & Done] j2004) also find that LMC X-l does 
not appear to follow the relation. However the range in disc tem- 
peratures and luminosities in their study, like the observations pre- 
sented here, is comparatively small. XTE 11650-500 is also in- 
cluded in their study, and although they find departures from the 
i'Disc oc T 4 law, these observations have large error bars and so 
are consistent still with it. 

The theoretical L — T relation depends on the mass of the 
BHXRB - which can reasonably be assumed to be constant dur- 
ing the outburst - the / co i, the colour temperature correction factor, 
the spin and also the inner radius of the disc, R[ n . The behaviour 
of the disc temperature and luminosity seen in Appendix Fig I A. II 
for the majority of the observations of the BHXRBs, indicates that 
over a large range in disc temperature and luminosity, both the f co \ 
and Ri n are relatively constant, at least in the most disc domi- 
nated states. There is some scatter around the best fitting relation 
which on the whole appears random with no clear secondary trend 
(see Fig. QJ. However, in the "spurs" at the beginning and ends of 
the outbursts, the deviation from the L — T relation is large, and 
we now investigate whether coherent variations either Ri n or / co i 
could cause them. 



6 DISC RADIUS 

Only when the inner radius of the disc, R[ n , is constant will the disc 
luminosity and temperature follow the expected L oc T 4 relation 
assuming that / co i is also constant. As we have shown in Section|5] 
that for the majority of the observations the expected L — T relation 



6 We note that t he not all of the b est fitting models appear appropriate 
for LMC X-l in IDunn et alj feOlOh . as all the observations have similar 
X-ray colours, yet some have no disc component a s seen in the Disc F rac- 
tion Lumin o sity D iagram (see lKording et al. 12006); Dunn et al. 1 20081) and 
IDunn et al] feoiol) for more details). This is likely to be the result of the 
similarity between a broken powerlaw model and a disc + powerlaw model 
when the disc does not dominate the spectrum. This is exacerbated in LMC 
X- 1 as the hard X-rays are also faint as a result of the large distance to the 
BHXRB, which make fitting the powerlaw difficult. 



Table 3. Disc Temperature Fits 



Object 


Exponent 


Points 


4U 1543-47 


5.76 ±0.76 


26 


4U 1630-47 


5.67 ±0.23 


245 


4U 1957+115 


4.10 ±0.24 


16 


GRO J1655-40 


4.16 ±0.14 


209 


GRS 1737-31 


- 




GRS 1739-278 


-0.59 ±0.02 


5 


GRS 1758-258 


4.80 ± 0.42 


2 


GS 1354-644 


- 




GS 2023+338 


- 




GX 339-4 


7.08 ±0.37 


197 


H 1743-322 


3.70 ±0.14 


99 


XTE Jill 8+480 


- 




XTE J1550-564 


4.26 ± 0.09 


109 


XTEJ1650-500 


24.64 ± 7.79 


28 


XTEJ1720-318 


5.21 ±0.19 


29 


XTE J1748-288 


3.19 ±0.16 


5 


XTE J1755-324 






XTEJ1817-330 


4.68 ±0.23 


88 


XTE J 1859+226 


4.11 ±0.34 


52 


XTE J2012+381 


4.28 ±0.55 


15 


LMC X-l 


-1.88 ±0.88 


14 


LMC X-3 


4.64 ±0.11 


109 


SAX 1711.6-3808 






SAX 1819.3-2525 






SLX 1746-331 


4.88 ± 0.40 


12 



Note 



The number of points are those used in the fitting, so with a Disc Frac- 
tion > 0.8. Notes: 1 - GRS 1739-278 has very few points for the fitting. 
2 - XTE J 1650-500 has a strange distribution of observations along the 
outburst, causing the erroneous fit. 3 - LMC X-l: small range in disc tem- 
peratures and luminosities and the slope is not well defined. 



is a good description of the behaviour, then we expect that the disc 
radius remains constant for the most disc dominated observations. 

To calculate the disc radius, we use the normalisation of the 
DISKBB model, JV, from XSPEC as this explicity includes the inner 
radius of the disc, R[ n . 



Rin 



D 



10 kpc 



cos 9, 



(2) 



where Diokpc is the distance in units if 10 kpc and 9 is the in- 
clination of the system. Where the inclination of the system is not 
known (Table Q3 we use 9 = 60° . It is possible that the discs in 
some of these BHXRBs are misaligned with resp ect to the binary's 
inclination (Maccarone 2002; Fragos et al. 2010). In Appendix Fig 



I A. II we show the inner radius against the disc temperature for each 
BHXRB. We do not include the effect of the uncertainties in the 
distance and inclination in the error bars shown. We also show in 
Appendix Fig lA. 11 a schematic diagram demonstrating the motion 
of the BHXRB through the R[ n — T plane as the outburst pro- 
gresses. As, on the whole, the disc temperature decays during the 
most disc dominated stages of the outburst, the track of the BHXRB 
through the _Ri n — T plane should be clear. 

As was expected from the behaviour of the disc luminosity 
and te mperature, and from the study of GX 339-4 bv lDunn et al.1 
(2008), the majority of points are at a relatively constant inner disc 
radius. We have chosen to plot the disc radius in kilometres on the 
primary x-axis rather than R as fractions of the gravitational ra- 
dius, _Rg, as in many of the BHXRBs the masses are not accurately 
known. The secondary x-axis shows the radii as a fraction of i? g 
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Table 4. Disc Inner Radius Fits 



Object 



Ri„(km) 



Rin/Rg 



4U 1543-47 
4U 1630-47 
4U 1957+115 
GRO J1655-40 
GRS 1739-278 
GRS 1758-258 
GX 339-4 
H 1743-322 
XTEJ1550-564 
XTE J 1650-500 
XTE J 1720-3 18 
XTE J1748-288 
XTEJ1817-330 
XTE J 1859+226 
XTEJ2012+381 
LMC X-l 
LMC X-3 
SLX 1746-331 



53.6 
20.9 
2.43 
21.8 
29.3 
4.90 
42.0 
22.4 
49.4 
18.9 
69.5 
28.2 
66.2 
35.5 
27.2 
38.2 
36.1 
6.7 



10.4 
3.6~ 
+0.07 
+2.1~ 
+0.6~ 
0.01 
+8.6~ 
+2.2 
+5.4~ 
+4.3~ 
16.7 
1.4 
+ 12.9 
+6.5~ 

1.4 
+5.6 
+2.l" 
+0.7 



8.7 
3.1 
0.07 
1.9 
0.6 
0.01 
7.1 
2.0 
4.8 
3.5 
13.4 
1.3 
10.8 
-5.5 
-1.4 
-4.9 
-2.0 
0.7 



3.86 
1.42 
0.16 
2.11 
1.98 
0.33 
4.90 
1.52 
3.15 
2.14 
4.70 

I. 91 

II. 1 
2.40^ 
1.84 
2.59 4 
2.44 
0.45 4 



+0.75 



1 


— 


02 


,+0. 


24 






-0 


21 


;+o 


.005 






-0 


.005 


+0 


21 






-0 


19 


;+° 


.04 






-0 


04 


;+° 


0009 






-0 


.0009 


1+ 1 


.00 






-0 


.83 


,+0 


15 






-0 


14 


!+° 


.34 






-0 


31 


+0 


.49 






-0 


40 


1+ 1 


.13 






-0 


91 


+0 


.09 






-0 


09 


+ 2. 18 








1.82 



+0.10 



+0.14 



The inner radius estimates do not include the uncertainties on the values 
or estimates on the distance or mass of the BHXRB (see Equation^). As a 
reminder, the inner radius of the disc is 6i? g for a non-rotating black hole, 
and li? g for a maximally rotating one. 

for comparison. We show in Table [2] the best fit inner radii for the 
observations with a Disc Fraction > 0.8. 

In some BHXRBs (4U 1957+115, GRS 1758-258 and 
SLX 1746-331) the disc radii are very small, less than 10 km or 
below 1-Rg. The innermost stable circular orbit (ISCO) for a black 
hole is 6-Rg for a non-rotating black hole, where R g = GM/c 2 = 
is the Schwarzschild radius. For a maximally rotating black hole, 
the ISCO can go down to 7? g . In these three BHXRBs the distances, 
masses and inclinations are not known, and so these effect the es- 
timates on the inner disc radius. From the normalisation, A/, the 
distance is directly proportional to the R[ n , the mass inversely pro- 
portional (when measured in units if R g ) and the inclination has a 
(cos dependence. The most change would arise if the incli- 

nation would increase, though increasing the distance or decreas- 
ing the mass would also increase the inner disc radius. However, 
not knowing the true values of the inner disc radius, it is difficult 
to determine which of these parameters should change and by how 
much. Therefore, although the estimates on the inner disc radius 
are smaller than physically sensible, it is likely to be the result of 
our incomplete knowledge of the system parameters. 

The "spurs" which were mentioned in Section[5]are also seen 
in the R[ n — T plots. If taken at face value, then they indicate that 
at the end of the outburst the inner disc radius decreases as the disc 
temperature rises, and vice versa at the beginning of the outburst. 
This behaviour does not appear to be physically meaningful, as the 
minimum radius measured for some of the observations fall well 
within the smallest stable orbit for a IOMq black hole. We dis- 
cuss these spurs and their possible causes further in the following 
section. 



7 THE SPURS - AN F CO l CONNECTION? 

As has been alluded to in the above sections, apart from the theo- 
retically expected L — T relation at constant disc radius 7?i n , the 



other main trend is perpendicular to the I/Disc or T 4 relation, lead- 
ing to an apparent decrease in the disc radius at the beginnings and 
ends of the outburst. These spurs a re seen in most of the BHXRBs, 
and w ere also seen in GX 339-4 in lpunn et alj I2008T) , XTE J1650- 
500 inlGierlinski & Pond J2004l) and GRO J1655-40 in lDone et all 
d2007h . 



The plots of Ri n — T in Appendix Fig. lA. ll show that the disc 
radius decreases to very small values at the beginnings and ends of 
the outbursts. It is the extent of the decreases, down to values well 
below 1-Rg, that lead us to investigate whether changes in / co i could 
be responsible. Although changing the system parameters can shift 
the location of the observations in the Ri n — T plane, they are un- 
likely to move the observations with the smallest calculated Ri n 
sufficiently far. We note that there are also points which occur at 
quite a distance from the main L — T relation, and also have very 
low disc fractions (< 0.5) and also have large uncertainties in the 
disc temperature and luminosity. These are unlikely to be explained 
by a variation in / co i, and are discussed further in Section[9] 

The limited low energy response of the RXTE PCA instru- 
ment may restrict the accurate fitting of a disc components when 
it does not dramatically dominate over the remainder of the con- 
tinuum. The calibrated range of the PCA starts at around 3keV 
whereas the disc temperatures peak at around 1 — 2 keV. As we 
are therefore fitting only one side of the disc component, as the 
temperature and the luminosity fade, the slight curvature could be 
difficult to accurately fit especially in short observations or ones 
which have a small number of signal counts. Although some of the 
spurs could arise from weaknesses in our data an alysis proceedure, 
as they have also been s een in other studies (e.g. Gierlinski & Pond 
l2004l : 1Pone et alj|2007t) it is l ikely that these weakness es are not 
the full explanation. However [Gierlinski & Pond J2004h only use 
the most disc dominated observations in their analysis. We have 
emulated this approach in this study when fitting lines or finding 
averages by selecting those observations with a very high disc frac- 
tion (> 0.8). However, the plots show all the observations which 
have a detected disc component. Therefore, our plots show these 
spurs, which may not have been shown in other studies where only 
the selected observations were plotted. 

W e also note that in a study of LMC X-3 by ISteiner et al.1 
(2010), the observations from RXTE PCA PCU-2 gave a very con- 
sistent value for the inner radius, a value which was also consistent 
with those obtained from other detectors (e.g. Suzaku, Swift and 
XMM-Newton). In our study we also find no evidence for the spurs 
seen in other BHXRBs in LMC X-3, and so their study does not 
help in clarifying whether the spurs come from instrumental effects. 
However the close correspondence of the R[ n between RXTE and 
the other detectors does indicate that in the disc dominated state, the 
spectral coverage of the PCA is sufficient to be able to accurately 
determine the disc parameters. 

In the theoretical L — T expectation, we have used a constant 
value for the colour correction, f co \ = 1.8 JShimura & Ta kahara 
Il995h . The departures from the L oc T 4 law could be the result 
of the variation of the value of / co i. The colour correction factor 
accounts for the change in the dominant emissi on process in the in- 
ner di sc, and therefore is a function of I/Disc- Gierlinski & Pond 
d2004l) showed that the / co i is approximately constant through- 
out the outburst (see als o lShimura & T akahara 1995; Merlo nTet al.l 
2000; IPavis et alT2 006). However these were only for the most disc 
dominated observations. Therefore it is possible that / co i variations 
could have occurred in the BHXRBs presented in these studies, but 
not be shown in the figures. If / co i is constant during the most disc 



© 0000 RAS, MNRAS 000, 000-000 



8 Dunn, Fender, Kording, Belloni & Merloni 



dominated phases of the outburst, and only those phases are shown, 
then any variation would not be detected. 

To account for these spurs, / co i would decrease at the begin- 
ning of the outburst, and then increase on the exit of the outburst as 
the BHXRB goes through the intermediate states (see the schematic 
in Appendix Fig, IA.lt . The motion off the T 4 relation is approxi- 
mately perpendicular. Therefore a simple change in f co \ would ex- 
plain the deviations, without needing any further variation (of inner 
disc radius, for example). It is of course possible that the disc ra- 
dius is not constant at the very beginnings and ends of the outburst. 
However, if the / co i is not constant then it will be difficult to deter- 
mine what the true R[ n is in these non dominant discs observed by 
RXTE. 

We show on the L — T plane in Appendix Fig. I A. 11 the the- 
oretically expected relation for / co i = 1.8 but also for a range of 
values for / co i between 1.6 and 2.6. The lower limit arises from 
the initial investigation into / co i by Shimura & Taka haral J 1995h . 
whereas the upper comes from the best characterised BHXRB in 
iGierlihski & Pond {2004). In cases where the values of the distance 
and mass used are such so that the theoretical L — T relation is a 
good match to the observed L — T relation, then the spurs, should 
they be present, mostly fall within this 1.6 < / co i < 2.6 range. In 
Appendix Fig. I A. 11 we also show a schematic which indicates the 
route taken by a BHXRB in this diagram. 

If these spurs are purely the result of changes in / co i, we can 
calculate the change in / co i required, <5/ co i, for the spurs to be part 
of the expected L — Rf n T 4 relation, under the assumption that the 
disc radius is constant. We assume that / co i = 1.8 when the T 4 
relation is followed, and so adjust the normalisation of Equation 
[T]so that the expected relation falls under the observed points at 
kT = 1 kcV. As this normalisation is affected by the distance, 
mass, and inclination, which in many BHXRBs are only estimates, 
this simplifies our approach, without affecting our conclusions on 
the variation of / co i . We also note that the spin of a black hole can 
affect the normalisation. We are currently assuming that the black 
hole is not rotating, but discuss spinning black holes in see Section 
[8] In Appendix Fig. lA. ll we show for each BHXRB the excess / co i 
required for the observation to lie on the T 4 relation, 5/ co i, against 
the powerlaw and disc fractions of the observation. In most cases 
this centred on Sf co \ = 0.0, which is by design, though where the 
L — T relation slope is very different from 4 (Section[5]l, then the 
position along the a>axis can vary. 

There appear to be three regions in the diagrams. The obser- 
vations with the largest disc fractions cluster around Sf co \ = 0.0, 
as defined by the normalisation adjustment mentioned above. The 
Sf coi remains almost constant at zero over around an order of mag- 
nitude change in the disc fraction. These are the observations which 
scatter around the theoretically expected T 4 relation, as the / co i is 
constant. 

As the disc fraction decreases the trend is for the observa- 
tions to move gradually towards progressively higher values of / co i, 
5/coi increases. These are the beginnings of the spurs, but are also 
visible as lopsidedness in the scatter around the most disc domi- 
nated observations in the plots of Luisc versus T. 

At around <5/ co i = 0.5 the trend in the observations flattens 
off, as the disc fraction approaches zero, resulting in large changes 
in / co i over small changes in the disc fraction. These observations 
are the ones from the spurs and extend up to 5/ co i ~ 1.0. The value 
of /coi required for these observations to lie on the T 4 relation be- 
come larger with very little change in the disc fraction. The x-axes 
of the figures has been truncated as the observations with very low 
disc fractions which lie well below the main cluster of points have 



up to S /coi ~ 20 (see Section [9}. These observations are unlikely 
to be explained by a varying / co i and hence we do not show them 

in the figure. 

The most recent investigation into / co i by iDone & Davis! 

(2008) shows that there is a positive correlation between the / co i 
and the mass accretion rate. The effect is stronger for a proportional 
counter array (e.g. the RXTE PCA) than for a Charge-Coupled De- 
vice (CCD) and shows that for an alpha disc, with a = 0.1, / co i 
can reach values of around 2.1 for accretion rates of 10 19 g/s (for 
a KERRBB disc model). This is lower than the / co i increase in- 
ferred in Appendix Fig. I A. II but links the accretion rate to the 
/coi- Although the e volution of the / co i over time indicated by 
IDone & Davis! d2008h is different, the clear link between the de- 
viations from L — T 4 and / co i suggests a link to the accretion 
rate. However, it must be noted that the color correction fraction 
and, in general, the observed properties of the high-energy tail of 
the disc emission are quite sensitive to the verti cal structure of the 
disc . In particular, as discu ssed in more detail in lDavis et al.1 J2005h 
and IDone & Davisl 1200 8D . the vertical dissipation profile may be 
very different in spectral states where a non-thermal (power-law) 
component is significantly detected, as a larger fraction of the total 
accretion power has to be released near or above the disc surface, 
leading to a possible increase in the estimated color correction fac- 
tors. 

Both a constant / co i and a constant Ri n are observed when 
the disc fraction is high (> 0.8), i.e. in those observations where 
the disc emission dominates over the powerlaw emission, and the 
disc parameters have been well determined. However, under the 
assumption of a constant / co i, as the disc fraction reduces the in- 
ferred Ri n decreases, which has a knock-on effect on the behaviour 
of the disc temperature with the luminosity. However, the drastic 
nature of the decrease in the Ri n is such, that a increase in / co i 
may be a more reasonable explanation. If the true underlying be- 
haviour of the BHXRB is that the inner radius remains constant, 
then an increase of the / co i could account for the majority of the 
spurs observed. 



8 DISTANCES, MASSES AND SPINS FROM DISC 
PARAMETERS 

During the most disc dominated parts of the outburst, the inner ra- 
dius of the disc and the / co i are approximately constant. The es- 
timated disc parameters (temperature and luminosity) depend on 
the physical parameters of the BHXRB system - the distance, mass 
and spin of the black hole. However, these parameters are constant 
for a particular BHXRB. Therefore, the shape that the observations 
make in the L — T plane are fixed, but their location within the 
plane could vary, depending on these parameters. 

Using the theoretical relation (EquationQJ and the data we are 
able to place limits on the distances and masses of the BHXRBs, 
in the case of a Schwarzschild black hole. From Equation Q] the 
theoretically expected L — T relation is of the form of 

^ifi = AM t\ (3) 

where M is the mass of the black hole. When fitting the trend in 
the L — T plane for the most disc dominated states, the form is (for 
constant R ln and / co i) 

^=B^-T C , (4) 
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where T is the observed temperature, B is derived from the normal- 
isation of the DISKBB component in XSPEC and D is the distance of 
the BHXRB. The dependence on M and D in these two equations 
are different. As these two Equations should be equal, assuming the 
fitted value of C = 4, then 



D 2 

AM = B±L. 

M 



(5) 



However, the slope of the fits to the L — T relation do not always 
end up with C ~ 4. Therefore, in order to remove this dependence 
we calculate the match at kT — 1 keV, which is close to the tem- 
peratures of the observed discs. 

Therefore, under the assumption that R[ n and / co i are con- 
stant, we can determine the distances and masses which are re- 
quired in order that the location of the expected L — T relation 
matches those which are observed. In cases where limits have been 
placed on either the distance or mass, then we are able to constrain 
the acceptable values for the mass or distance respectively. These 
loci of points in the distance-mass diagram are shown in the Ap- 
pendix, Fig. I A. 1 1 by the blue line, and we also show the current 
best estimates on the distances and masses and their uncertainties, 
where they exist. In some cases (e.g. GRS J1739-278, XTE 11650- 
500 LMC X-l) the fitted slope is very different to 4, and so there in 
these cases, these loci are not reliable. 

In some cases the current best observational estimates on the 
distances and masses do coincide with the estimates from this 
work (e.g.4U 1543-47, GX 339-4, XTE J1720-318). However in 
many cases there is no overlap between the observed estimates on 
the distance and mass and those calculated here (e.g. GRO 1655- 
40, XTE J 1550-564). Although at face value, cases where there 
is no overlap would allow the distance and mass estimates to be 
refined, it is not quite that simple. These constraints are for a 
non-rotating black hole, and there is significant evid ence that at 
least s ome black holes have si g nificant spin (see e.g. Miller et all 
20091: iMcClintocketalJ [2006; Midd leton et id] l200d " and also 
Fender et al. 2010 and references therein). 



The normalisation of the theoretically expected relation be- 
tween LDisc and T changes w hen the black hole is m aximally 
spinning (a* = 0.998) (see e.g. lGierlinski & Dondl20q4). We use 
a very simple parameterisation from Ma kishima et al. who 
include first order effects of the black hole spin on the theoretically 
expected relation in terms of 

a = R in /3Rs, 

the ratio of the inner disc radius to the Schwarzschild radius (a = 
1 and a — 1/6 for a non rotating and maximally rotating black 
hole respectively). This appears as an a 2 term in their version of 
Equation[T] We add this correction factor into EquationQ] and also 
take into ac count the changes i n the general relativistic correction 
factors from lZhang et al] J 19971) for rotating black holes. 

Therefore, we show in Appendix, Fig lA.ll as well as the range 
of distances and masses allowed for a non-rotating black hole, we 
also show those for a maximally rotating black hole (a* = 0.998, 
red line) and one for a* = 0.5 (green line). This results in an area 
in the D — Al plane in which both values are allowed. We note, that, 
a more accurate investigation would start with a more appropriate 
model for the disc emission including the relativistic effects of the 
black hole spin (e.g. KERRBB). We have, as yet, not re-fitted all 
our results with such a model. These plots show that many of the 
BHXRBs whose distances and masses did not match the loci for a 
non-rotating black hole, do match if the black hole is rotating (e.g. 
GRO J1655-40, XTE J1550-564, XTE J1650-500, LMC X-3). The 




1.2 1.4 1.6 

Log(Break En) (keV) 

Figure 2. The distribution of the break energy for the broken power- 
law fits using two different binning levels. The gap between 1.4 < 
log 10 i?Break < 1-6 arises from the cross over of the PCA and HEXTE 
instruments. 



eight BHXRBs which have both mass and distance estimates, are 
all consistent with spin values < a* < 0.998, or, alternatively the 
inner radii are consistent with 7? g < R[ n < 6R S . In fact, none of 
the estimates on the spin of the black hole require high spin values 
(a ~ 1). This also indicates that there are no counter-rotating discs 
in these BHXRBs. As the estimated spins fall in the range expected 
for Kerr black holes, this indicates that the position of 7?i n in the 
disc dominated states is mainly determined by strong gravitational 
effects. 

For many of the BHXRBs in this sample, either the mass or the 
distance or both are unknown. If only one is known, then limits on 
the other can be placed from the range allowed by the black holes 
spin. However, in many cases these limits are not very constrain- 
ing (e.g. GRS J1739-278, XTE J1720-318). If both are unknown, 
then the best that can be obtained is a lower limit on the distance 
(assuming a reasonable lower bound for the black hole mass). 



9 BROKEN POWERLAW - DISC MODEL DEGENERACY 

The limitations of the spectral response of the RXTE PCA make the 
unequivocal detection of a non-dominant disc very difficult. When 
the disc's spectral component begins to rise above the powerlaw, 
the difference in the xt between the broken powerlaw and disc + 
powerlaw model is very small. It is therefore difficult to say ab 
initio which of the two models is the m ost appropriate to select 
for the observation (see lDunn et al.ll2010l) . In our selection proce- 
dure we select merely on the lowest reduced x 2 , with some restric- 
tions on the model parameters. We now investigate how adapting 
the model selection proceedure affects the numbers and parameters 
of the discs detected. 

As the broken powerlaw model could mimic a disc model, 
especially when the break energy is low (< lOkeV). We show 
in Fig. [2] the distribution of the break energies in the observations 
best fit by the broken powerlaw model, using two different binning 
schemes. The gap around log 10 .Ebreak = 1.4 is the result of the 
crossover between the PCA and the HEXTE instruments at 25 keV. 

The main cluster of points occur at low energies, with another 
smaller cluster at around log 10 inbreak = 1.8. These observations 
with high break energies (-EWcak ~ 60 keV) are likely to be re- 
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Figure 3. The variation of the disc luminosity with the disc temperature 
for GX 339-4 and H 1743-322 when penalising the \ 2 of broken powerlaw 
models which have a break energy < 15keV. For comparison plots see 
Appendix Fig. lA. II 

liable fits to a true break in the spectrum. However, those which 
fall below around £Wcak ~ 10 keV could be the result of a non- 
dominant disc mimicking a broken powerlaw 

We redo the model selection proceedure but penalise the bro- 
ken powerlaw models where the break energy is lower than 15 keV. 
In some cases only one of the two powerlaw fits (with and without 
line) has a break energy < 15 keV, and then the remaining bro- 
ken powerlaw fit is still allowed. However, this only occurs in rare 
cases. The other two remaining options are then single powerlaw or 
disc+powerlaw models (each with and without a line). We expect 
that this will have most effect in the intermediate states, where non- 
dominant discs are expected, rather than in the most disc dominated 
states where L oc T 4 . 

In Fig.[3]we show the L - T plane for GX 339-4 and H 1743- 
322 as an example to show the increase in the number of observa- 
tions which have a low disc fraction. These extra observations ap- 
pear below the T 4 relation, and the spurs appear to merge in with 
them. Most of these observations have large uncertainties in both 
the disc temperature and luminosity. 

Doing the converse on the selection procedure - favouring the 
broken powerlaw models in these situations - removes most of the 
observations with a disc detection at very low disc fractions (lower 
right in the L — T plane) in Appendix Fig. IA.1I Although this is 



may be a more conservative selection procedure (no discs detected 
when they may be uncertain), it is less informative on the behaviour 
of the BHXRBs on the transitions between the different states. 

It is therefore difficult to determine what the appropriate se- 
lection procedure is when broken powerlaw and disc models are 
both good fits to the data. The spectral resolution and low energy 
range of the RXTE PCA are very limiting in this case. Future inves- 
tigations or instruments may aid in pin-pointing the disc behaviour 
in the transition regions 



10 POWERLAW EVOLUTION 

Using the Disc Fraction Luminosity Diagram (DFLD) to investi- 
gate the behaviour of a BHXRB during an outburst restricts what 
information can be extracted about the variation of the powerlaw 
component, especially in the hard/powerlaw-dominated state. In 
this state, the variation in the powerlaw slope or break energy do 
not effect the disc fraction, and so all observations fall on a single 
line. 

In Appendix Fig. |A.2| we show the DFLDs for the X-ray bi- 
nary in question, where the colourscale shows the variation of the 
powerlaw slope (below the break if it is a broken powerlaw). We 
show in the neighbouring panel, the powerlaw slope against the to- 
tal luminosity for those observations with a disc fraction of < 0.2 
(powerlaw fraction > 0.8). This allows the change in the powerlaw 
slope to be tracked in the powerlaw dominated state. In a number of 
binaries there are insufficient observations in the hard state to de- 
termine any trend with time. Also, 4U 1630-47, the variation of the 
powerlaw slope appears complex, with no easily discernable global 
trend. However, the outburst structure in this BHXRB is als o com- 
plex and so this variation is expected (see lDunn et alj 2010). 

In the majority of BHXRBs, in the low luminosity "stalk", 
the powerlaw slope increases as the luminosity falls - the spectrum 
softens. This has been seen in the HIDs of the BHXRBs before, 
as a change in the X-ray colour. The re-emergence of the disc at 
very low luminosities has been observed in deep pointed observa- 
tions and may also pl ay a role in the softe ning of the spectrum at 
low luminosities (see Cabanac et al. 2009). However, at these low 
luminosities the effects of the Galactic Ridge Emission (GRE) play 
a role. None of the BHXRBs were fitted with a model which takes 
into a ccount the effects of the GRE in their vicinity l lDunn et all 
2010). At low luminosities, the GRE can have an appreciable effect 
on the shape and flux of the spectrum. Howev er, the curvature i n 
the stalk was seen in the study of GX 339-4 bv lDunn etall d2008l) . 
where the GRE was added to the model spectrum as a fixed com- 
ponent. Therefore only part of the softening at low luminosities can 
be explained by the GRE. 

However, at the top of the powerlaw dominated state in some 
BHXRBs, there is an increase in the powerlaw slope. The increase 
in r has been observed as the BHXRB enters the soft or disc- 
dominated state. However, in these Figures, the V increases far be- 
yond what has been ob served in other st udies of these BHXRBs, us- 
ing the same data (e.g jMotta et alj2009h . The likeliest explanation 
is that the broken powerlaw is accounting for a rising disc, which 
is not being well fit by a disc model. What can also be seen is that 
these softer powerlaw slopes are from broken powerla ws, and they 
have a comparatively low break energy. As noted in iDunn et al.l 
(2010) and Section|9]there is a possibility for the broken powerlaw 
to mimic the disc (and powerlaw) model. It is probable that as the 
disc rises in luminosity the limited spectral range of the RXTE PCA 
means that the curvature of the disc cannot be determined, and the 



© 0000 RAS, MNRAS 000, 000-000 



A Global Study of the Behaviour of Black Hole X-ray Binary Discs 1 1 



broken powerlaw resulted being a better fit. Restricting the power- 
law break to being above the peak mentioned in Section [9] would 
prevent this occurring. However, the accuracy of the fitted disc pa- 
rameters is not clear. Therefore, the softening of the powerlaw on 
the transitions to the disc-dominated state is likely to be the result 
of the limitations of the RXTE PCA. 



11 SUMMARY 

We have investigated the behaviour of th e disc and powerla w com- 
ponents in the 25 BHXRBs presented in lDunn et alj d2010h . In the 
majority of BHRXBs in which at least most of an outburst has been 
observed, the disc luminosity scales close to T 4 in the most disc 
dominated observations. This behav iour had been seen in other 
studies (e.g. lGierliriski & Pond 120041) . The scaling of T 4 implies 
that both the disc's inner radius, R ln and the colour correction fac- 
tor, /coi are relatively constant in the most disc dominated states. 
A number of BHXRBs do not show a clear T 4 relation, but these 
could be the result of the limitations of the model fitting routine or 
the frequency of observations. 

However, in observations where the disc is no longer overly 
dominant, there are deviations from the T 4 law. If interpreted as 
changes in the disc's inner radius, these deviations imply that the 
disc is moving inwards at the end of an outburst, and outwards at the 
beginning of the outburst. Although we do not rule this behaviour 
out, it seems an unlikely scenario. If these deviations are attributed 
to changes in the colour correction factor, then / co i rises as the disc 
fraction decreases. 

There are a number of observations in which the disc param- 
eters determined are unlikely to be explained by reasonable values 
for the disc radius or the / co i. The spectral fits for these observa- 
tions tend to have \ values which are very similar to those for the 
broken-powerlaw model, which makes selecting the most appropri- 
ate model difficult. This also makes determining the true behaviour 
of the disc temperature, radius and / co i in these intermediate states 
difficult. 

The luminosity of the disc in Eddington units can be calcu- 
lated from theoretical arguments from the BHXRB parameters and 
the disc temperature. When calculating the observed luminosity the 
distance of the BHXRB system also enters the calculation. The 
combination of these two calculations allows the ratio of D/M to 
be estimated from the L — T relation of the BHXRB disc. We have 
therefore placed limits on the range values of D and M values al- 
lowed for different values of the spin of the black hole. 
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APPENDIX 

For each BHXRB in our sample we show in Fig. lA.ll TOP LEFT the 
disc temperature as a function of the disc luminosity, along with 
the theoretically expected relation for the case that the BH is not 
rotating for a variety of values for / co i. The f co \ = 1.8 line is 
depicted thicker than the rest, which increment in 0.2 intervals. The 
line which best fits the most disc dominated spectra is shown, as 
well as its slope. We show TOP RIGHT the region in the distance- 
mass plane which is allowed by the theoretical expectation of the 
luminosity- temperature relation and the observed data points. The 
best determined values for the distance and mass along with their 
uncertainties are also shown where available. 

In the BOTTOM LEFT we show the variation of the disc radius 
(Rin) with the disc temperature. The colourscale is the disc frac- 
tion. The line is the average of the disc radii for the observations 
which have a disc fraction > 0.8. The variation of the / co i with the 
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Figure A.l. Schematic figures of LEFT: Disc luminosity against disc temperature, MIDDLE: Disc radius against disc temperature and RIGHT: the excess (5/ co i 
required against the disc fraction. The arrows show the motion through the diagram. In the right-hand figure, the offset between the inward and outward tracks 
are for clarity only. 




Figure A.l. (cont) 4U 1543-47 TOP LEFT We show the Luminosity of the disc as a function of the temperature of the DISCBB component in XSPEC.The colour 
scale is the disc fraction of the observation. The theoretical L — T relations are shown by the sets of dotted line for a number of values of / co i . The solid line 
shows the fit to the most disc dominated states (Disc Fraction > 0.8). TOP RIGHT We show the range of distances and masses allowed, if the observations are 
to match the theoretical relation. The blue line is for a non-rotating black hole, and the direction for increasing a is shown. The dashed and dotted lines show 
the current values and uncertainties on the masses and distances where available. BOTTOM LEFT We show the variation of the inner disc radius (Ri n ) with 
the disc temperature. The colour scale is the disc fraction of the observation. The solid line shows the average Ri a of the most disc dominated states (Disc 
Fraction > 0.8). BOTTOM RIGHT The excess the colour correction required, <5/ co i, with the disc fraction. The disc temperature is the colour scale. 



disc fraction is show in the BOTTOM RIGHT. The disc temperature 
is the colour scale. 

By comparing all three of the scatter plot figures, the be- 
haviour of the disc in the BHXRB in the outburst becomes clearer. 
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Figure A.l. (cont) 4U 1630-47 
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Figure A.l. (cont) 4U 1957+115 
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Figure A.l. (cont) GRO J1655-40 
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Figure A.l. (cont) GRS 1737-31 
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Figure A.l. (cont) GX 339-4 
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Figure A.l. (cont) H 1743-322 
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Figure A.l. (cont) XTE J1550-564 
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Figure A.l. (cont) XTE J1650-500 
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Figure A.l. (cont) XTE J1720-318 
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Figure A.l. (cont) XT J1748-288 
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Figure A.l. (cont) XTE J1817-330 
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Figure A.l. (cont) XTE J1859+226 




Figure A.l. (cont) XTE J2012+381 
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Figure A.l. (cont) LMC X-l 
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Figure A.l. (cont) LMC X-3 
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Figure A.l. (cont) SLX 1746-331 
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Eisjire^ A.2. The DFLDj^ith Jhe^colour scale showing the powerlaw slope. Also shown in the side panel is the powerlaw slope (below the break where 
a^ropna%^^b^ffl?^servation^w?8i a disc fraction < 0.2, along with the break energy for the broken powerlaw (where it occurs). TOP: 4U 1543-47, 
middle: 4U 1630-47, bottom: 4U 1957+115. 
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Figure A.2. (cent) TOP: GRO J1655-40, middle: GX 339-4, bottom: H 1743-32^ ^ ^ MNRAS 000, 000-000 
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( (cont) TOP: XTEJ1 1184480, middle: XTE J 1550-564, bottom: XTE J1650-500. 
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Figure A.2. (cont) TOP: XTE J1720-318, MIDDLE: XTE J1748-288, BOTTOM: XTE J18 
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Figure A.2. (cont) TOP: SAX 1711.6-3808, middle: SAX 1918.3-2525, bottom: SLX l^jgl,^ MNRAS 00Q 000 _o 00 
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